One-loop self-energy correction in a strong binding field 
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I. INTRODUCTION 

Calculations of the one-loop self-energy correction to all 
orders in the parameter Za (Z is the nuclear charge number 
and a is the fine-structure constant) have a long history. The 
first correct evaluation of this correction was performed for 
several high-Z ions by Desiderio and Johnson (jj using the 
method proposed by Brown, Langer, and Schaefer |2||. An- 
other, much more accurate and powerful method was devel- 
oped by Mohr 1 3], who carried out high-precision calculations 
of this correction in a wide range of Z for the ground and the 
first excited states of H-like ions fl, IH |(J . Various extensions 
of this method provided highly accurate evaluations of the 
self-energy correction for higher excited states 0,01, for the 
extended-nucleus Coulomb potential |9], and for very small 
nuclear charge numbers BlQt ll ll fl2ll . Indelicato and Mohr 
JT3I [fl presented an important modification of the method, 
in which renormalization is performed completely in coordi- 
nate space. 

A different method for evaluation of the self-energy cor- 
rection, which can be conventionally termed as the potential- 
expansion method, was introduced by Snyderman and Blun- 
dell 0E0. Various numerical schemes based on this 
method were presented by other groups fTsl fT^il . 

There are also other methods developed for evaluation of 
the self-energy correction which have been less widely used 
so far. A noncovariant method of the so-called partial-wave 
renormalization was developed by Persson, Lindgren, and Sa- 
lomonson 1 20] and by Quiney and Grant 12 U |22| |. Another 
method proposed by Labzowsky and Goidenko 12311 is based 
on the multiple commutator expansion of the general expres- 
sions. 

Closely related to the self-energy is the other dominant 
QED effect, the vacuum-polarization. The first evaluations of 
this correction to all orders in Za were performed by Soff and 
Mohr 1 24] and by Manakov, Nekipelov, and Fainstein 12511 . 
More accurate calculations of the vacuum-polarization correc- 
tion were carried out later by other groups j26ll27il . 

Evaluation of the self-energy correction for a tightly bound 
electron is nontrivial, to a large extent, due to the fact that this 
correction involves the Dirac-Coulomb Green function that is 
not presently known in the closed analytical form (contrary to 
the nonrelativistic Coulomb Green function). Consequently, 
the self-energy correction is expressed as an infinite expan- 



sion over the angular momentum of the virtual photon (or, 
equivalently, the total angular momentum of the virtual elec- 
tron states j = \k\ — 1/2, where k is the relativistic angular- 
momentum parameter of the Dirac equation). This expan- 
sion (further referred to as the partial-wave expansion) greatly 
complicates calculations of the self-energy corrections. 

In the method by Mohr [3], the summation of the partial- 
wave expansion was performed numerically before integra- 
tions over radial coordinates. A large number of terms in- 
cluded into the summation (~ 10 4 ) and usage of the quadruple 
arithmetics ensured a high accuracy of the numerical results 
obtained but made the computation rather time consumin g. I n 
the extension of this method by Jentschura et al. II (ill llll 211. 
several millions of expansion terms included into computa- 
tion were reported, which became possible due to an elaborate 
convergence-acceleration technique developed by the authors 
and an extensive usage of modern parallel computer systems. 

On the contrary, calculations based on the potential- 
expansion method fla, flTI. Il8l fl9t1 are usually performed with 
much smaller numbers of partial-wave expansion terms ac- 
tually included into the computation (~ 15 — 40). This is 
achieved (j) by employing a more complete set of renormal- 
ization terms that are calculated separately in a closed form, 
(z'z) by performing the radial integrations before the partial- 
wave summation (for the discussion of how this influences the 
convergence rate see Eqs. (1), (2) of Ref. 111411 and the related 
text there), and (Hi) by using extrapolation to estimate the con- 
tribution of the tail of the expansion. The price to pay is a 
more complex structure of the subtraction terms (especially, 
in coordinate space) and the necessity to keep the accuracy 
of numerical integrations well under control for each partial- 
wave term, in order to provide a reasonable extrapolation for 
the tail of the expansion. Still, the method is computationally 
very cheap and can be directly generalized for calculations of 
higher-order QED diagrams, where the self-energy loop en- 
ters as a subgraph. These advantages have determined the fact 
that most calculations of higher-order self-energy corrections 
have been performed by extensions of the potential-expansion 
method up to now. 

The one-loop self-energy correction is traditionally repre- 
sented in terms of the dimensionless function F(Za), which 
is connected to the energy shift (in units H = c = m = 1) by 



AE=«^F(Za), 



(1) 
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where n is the principal quantum number. Practical calcula- 
tions performed within the potential-expansion method in the 
Feynman gauge show that the general behavior of individual 
partial-wave expansion contributions to the function F(Za) 
roughly follows the dependence 



10(Za) 2 |«| E 



(2) 



This makes clear that, while the nominal rate of convergence 
of the partial-wave expansion is always close to |k| -3 in this 
method, the actual convergence is governed by the parame- 
ter n 3 /(Za) 2 , whose numerical value can be rather large for 
excited states and small nuclear-charge numbers. Taking into 
account that the extension of the partial-wave summation be- 
yond the typical limit of |«| = 30 — 40 leads to serious tech- 
nical problems within the numerical scheme employed, we 
conclude that the parameter n 3 /(Za) 2 defines the region of 
the practical applicability of the potential-expansion method. 

Similar situation persists in calculations of self-energy cor- 
rections to higher orders of perturbation theory. In such cal- 
culations, the convergence of the partial-wave expansion also 
worsens with decrease of Z and increase of n. In particu- 
lar, a slow convergence of this expansion turned out to be the 
factor limiting the accuracy in evaluations of the self-energy 
correction to the Is and 2s hyperfine splitting in low-Z ions 
l28l I29I. l3(ill . This convergence also posed serious problems 
in calculations of the self-energy correction to the bound- 
electron g factor in light H-like ions 

The convergence rate of the partial-wave expansion be- 
comes most crucial in the case of two-loop self-energy cor- 
rections, for which the summation should be performed over 
two independent expansion parameters, both of which are un- 
bound I35ll36ll . A calculation of the two-loop self-energy cor- 
rection for very low nuclear charge numbers (and, specifically, 
for hydrogen) is a challenging problem, which apparently can- 
not be solved within a straightforward generalization of the 
potential-expansion method. (The present status of calcula- 
tions of the two-loop self-energy correction can be found in 
Ref. I37IO . One of the problems to be solved to this end is 
to find a way to improve the convergence properties of the 
partial-wave expansion. 

The goal of the present investigation is to formulate a 
scheme for evaluation of the one-loop self-energy correction, 
which yield the fastest convergence of the partial-wave expan- 
sion among the methods reported so far in the literature. 



II. FORMALISM 

The energy shift of a bound electron due to the first-order 
self-energy correction is given by the real part of the expres- 
sion 

AE = 2iaJ duj J dxi dx 2 x i2 ) 

x^(xi)a M G(e a - w,xi,x 2 ) a v V>a(x 2 ) 
-5m /dx^(x)/3^ a (x), (3) 



where a M = (1, a), a and ft are the Dirac matrices, 
G(w,xi,x 3 ) = [u - H(l - i0)]-\ H = Ho + V(x), 
Ho = a ■ p + ft is the free Dirac Hamiltonian, V(x) is a 
local potential (not necessarily the Coulomb one), and 5m is 
the mass counterterm. D^ v is the photon propagator defined 
in the Feynman gauge as 



D""(w,x 12 ) = 9 



_ M „ exp(i\/^ 2 + i0x 12 ) 



Alt Xl 



(4) 



where x\% = |xia| = |xi — x 2 |, and the branch of the square 
root is fixed by the condition Im(V^ 2 + i0) > 0. In Eq. Q 
it is assumed that the unrenormalized part of the expression 
and the mass counterterm are regularized in a certain covariant 
way and that the limit removing the regularization is taken 
after the cancellation of the divergent terms. 

Ultraviolet divergencies in Eq. can be conveniently iso- 
lated by separating the first two terms in the expansion of the 
bound-electron propagator G in terms of the binding potential 
V, 

G(£,x x ,x 2 ) = G (0) (£;,x 1 ,x 2 ) + G (1) ( J E;,x 1 ,x 2 ) 

+ G( 2 +)(£,x 1 ,x 2 ), (5) 

where G (0) = [w - H (l - i0)] _1 is the free Dirac Green 
function, G^ is the first-order expansion term 

G (1) (£,xi,x 2 ) = J dzG (0) (£,xi,z)U(z)G (0) (£;,z,x 2 ), 

(6) 

and G' 2+ - ) is the remainder. The three terms in Eq. (0, after 
substitution into Eq. (|3}, lead to the separation of the self- 
energy correction into the zero-potential, one-potential, and 
many-potential parts |15]: 

AE = AE zero + AE onc + AE many , (7) 

with the mass-counterterm part naturally ascribed to the zero- 
potential term. Converting the first two terms into momentum 
space and cancelling the ultraviolet divergences, one obtains: 



AE ZCYO 



dp 



^ a (p)sg ) ( £a ,p)V„(p), (8) 



/" dp! dp 2 - 

xr^(£ a ,pi;e a ,p 2 )U(q) Va(P2), (9) 

where q = p x - p 2 , ip a {p) = V4(p)7°, and (p) and 
r^(pi,p 2 ) are the renormalized free self-energy and vertex 
functions (for their exact definition and calculational formulas 
see, e.g., Ref. ifTl l. 

The many -potential term is represented by the following ex- 
pression 

A.E many = 2ia J duj J efoci ebc 2 D^(uj, x i2 ) 

x Vl(xi) G (2+) (e a -w,xi,x 2 )a„ -0a(x 2 ), 

(10) 
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where G*( 2 +) = G - G (0) - G (1) and the contour G of the 
uj integration does not necessarily go along the real axis but 
can be chosen differently in order to simplify the numerical 
evaluation of this expression. In our approach, we employ the 
contour Clh that consists of the low-energy part (Cl) and the 
high-energy part (Ch) and is similar to the one introduced in 
our previous work IU9II . The low-energy part of the contour 
Cl extends from — Eq — iO to — iO on the lower bank of the 
branch cut of the photon propagator and from +i0 to Eq + iO 
on the upper bank of the cut. In order to avoid appearance of 
poles of the electron propagator near the integration contour, 
each part of Cl is bent into the complex plane if the calcula- 
tion is performed for excited states. (The analytical structure 
of the integrand and a possible choice of the contour are dis- 
cussed in Ref. 1 19].) The high-energy part of the contour is 
Ch = (eo — ioo. e — iO] + [eq + iO, Eq + ioo). The parameter 
Eo separating the low- and the high-energy part of the contour 
is chosen to be £o = Za £ a in this work. (It is assumed that 
the condition £ a — e 1s < e is fulfilled for the states under 
consideration, where £i s is the ground-state energy.) 

Due to a lack of a closed-form representation for the Dirac 



Coulomb Green function, the evaluation of the many -potential 
term has to be performed by expanding G (and, therefore, 
G( 2+ )) into eigenfunctions of the Dirac angular momentum 
with the eigenvalue k. As discussed in Introduction, the con- 
vergence rate of the resulting partial-wave expansion is of cru- 
cial importance for the numerical evaluation of the self-energy 
correction. 

Until this moment, our description closely followed the 
standard potential-expansion method |15]. We would like 
now to modify this method in order to achieve a better con- 
vergence of the partial-wave expansion in the many-potential 
term Ai? many . To this end, we look for an approximation 

G { a +) to the function G< 2+ > that fulfills the following require- 
ments: (f) it can be evaluated in a closed form (i.e., without the 
partial-wave expansion) and (if) the difference G ( ^ 2+ - ) — Gi 2+ ^ 
inserted into Eq. dlOt yields a rapidly converging partial-wave 
series. 

We start with the expansion of the bound-electron Green 
function in terms of the binding potential, 



G(£,x 1; x 2 ) = G (0) (£, Xl ,x 2 )+ J dzG {a \E,x 1 ,z)V(z)G {0 \E,z, X2 ) 

+ J dz x dz 2 G ( °) (E, X! , zi) V( Zl ) G<°> (E, z u z 2 ) V(z 2 ) G<°> {E, z 2 , x 2 ) 



(11) 



r 



It is well known that the dominant contribution to radial in- 
tegrals like those that appear in Eq. (II Oi originates from the 
region where the radial arguments are close to each other, 
xi w x 2 . This region is also responsible for the part of 
the partial-wave expansion of the Green function that has the 
slowest asymptotic convergence in 1/|k| |3J. In this region 
the commutators of the potential V with the free Green func- 
tion G^ ) are small and can be neglected, which corresponds 
to expanding V(z) in a Taylor series around z = xi (or x 2 ) 
and keeping only the first term. Commuting V out to the left 
in Eq. Jl li and repeatedly employing the identity 

J dzG {0) (E,x 1 ,z)G i ° ) (E,z,x 2 ) 

= -^G(°)( J B,x 1 ,x 2 ), (12) 

we obtain the approximation G a to the bound-electron Green 
function G, 

xi,x 2 ) 



-nzi)^G(°)(£, Xl) x 2 ) 

+ V 2 (x 1 )-^G^(E,x 1 , X2 ) + 



This expansion has a form of the Taylor series and can be 
formally summed up, yielding 



G„(£,x 1 ,x 2 ) = G< >(.E + n,xi,X2), 



(14) 



where O = —V(x\) = Za/x\. Commuting V out to the 
right in Eq. M \\ . we obtain the same representation for G a 
but with SI = Za/x 2 . 

It should be noted that the idea of commuting the poten- 
tial V outside in the one-potential term was first proposed 
by Mohr ||3j], who proved that this procedure does not in- 
fluence the asymptotic ultraviolet behavior of this term (we 
recall that ultraviolet divergences originate from the region 
X! w x 2 in configuration space). Later, it was also demon- 
strated 1 13, 14] that all ultraviolet divergences in the one-loop 
self-energy correction could be identified by isolating several 
first terms of the power-series expansion of the potential V 
and the reference-state wave functions ip a around the point 
xi = x 2 . 

Expression dl4> yields an approximation for the bound- 
electron Green function that has a form of the free Green func- 
tion with a shifted energy argument. Taking into account that 
the free Green function is known in a closed form [ 3 ] 







[(-- 


X 12 / 



i a ■ xi2 + (3 + E 



(13) 



cxp[-ca;i 2 ] 

47T219 



(15) 
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(c = yl — E 2 ), we can employ this expression for the eval- 
uation of G a . 

An analogous to Eq. H41 approximation for the function 
G ( - 2+ - ) is obtained by subtracting the first two terms of the Tay- 
lor expansion from G a , 

G a 2+) (E, Xl , x 2 ) = G(°) (E + n, xi, x 2 ) - G<°) (£, X! , x 2 ) 



O^G(°)(^ ) x 1 ,x 2 ) 



(16) 



According to the derivation, the functions G a (E, Xi,X2) 

(2+) 

and Go (i?,xi,X2) approximate, correspondingly, 
G(£',Xi,X2) and G^ 2+ - ) (i?,xi,X2) in the region where 
xi w x 2 . This means, in particular, that instead of the 
original expression for 57 in Eq. (I14> . SI = Za/xi, one can 
use its arbitrary symmetrization with respect to x\ and x 2 - 
In our actual calculations, the following choice of SI was 
employed 



n = 



2Zc 



Xl + X 2 



(17) 



which turned out to be more convenient from the numerical 
point of view. 

We now use the approximate expression for the Green func- 
tion obtained above in order to separate the many-potential 
term dlOl into two parts, one of which contains G a 2+ ^ in- 
stead of G^ 2+ ) and is evaluated in a closed form in configura- 
tion space, whereas the remainder is calculated by summing a 
rapidly-converging partial-wave series. Bearing in mind that 
the partial-wave expansion for the low-energy part of Eq. (II 01 
is already converging very fast (if the parameter £o of the in- 
tegration contour Clh is chosen as described above), we ap- 
ply this separation to the high-energy part only. The many- 
potential term is thus written as a sum of the subtraction and 
the remainder term, 



AE — AE suh 



AE, 



rcmd 



(18) 



The subtraction term is obtained from the high-energy part 
of Eq. (US) by the substitution G< 2+ ) -> G { a 2+) . Its explicit 



expression in the Feynman gauge reads 



A£; sub = i« 
many ^ 



a f a a exp(i|w|xi 2 ) |, > 
dw / dxi dx 2 Va( x i) 

C H J ^12 

Gi 2+ )(e a -w,xi,X2)a^^ a (x 2 ). (19) 



The remainder term is obtained from Eq. (1101 by applying the 
substitution G (2+) — ► G (2+) - G a 2+) in the high-energy part. 



Calculational formulas for the remainder term AE, 



many 



are 



obtained by obvious modifications of the corresponding ex- 
pressions for the many -potential term that can be found, e.g., 
in Ref. d. In order to obtain the subtraction term in a form 
suitable for the numerical evaluation, one has first to perform 
the angular part of integrations over Xi, X2 analytically. To 

(24-) 

do so, we utilize the fact that both Ga and the scalar part 
of the photon propagator depend on angular variables through 
X12 only. Their product can be written as 



G( 2 +)( £a -^x 1 ,x 2 ) 



exp(i \u\ xn) 



X12 

T\ i a • X12 + T2 + Tz ■ 



(20) 



Here, Ti = J 7 i(ui,xi ) x 2> ^) are scalar functions depending 
on the radial variables through x\, x 2 , and £ = xi • X2 only, 
where x = x/x. Explicit expressions for Ti are immediately 

obtained from the definition of G a 2+ ^ dl6> and the expression 
for the free Green function G^ -* J15I . The functions T% can 
be expanded over the set of spherical harmonics by 

J : i (w,xi,x 2 ,0 = 4n^V l ^\uJ,Xi,x 2 )Yi m (-k 1 )Y l l l (x. 2 ) , 



(21) 



where 



1 r 1 

-J d^^i{u, Xl ,X2,0Pi(0 (22) 



and is a Legendre polynomial. 



Substituting Eq. J2 1 b into Eq. d!9l > and performing simple angular-momentum algebraic manipulations, we obtain 



AE™* ny = 2ia duJ d Xl dX2 d£,{x 1 X2) 2 {T 1 {LO,X ll X2^)9a{xi)fa{x2)[xiP 1 {0-X2Pl a {0} 

Jc H Jo J-i 1 

+ J r l(uj,xi,x 2 ,0 fa{xi)ga(x 2 ) [x 2 Pj a (0 ~ x x Pi a (£)] + 2T 2 (w, xt , x 2 , f) [g a {x l ) g a {x 2 ) Pi a (0 

-fa(xi)f a (x 2 )P Ja (0] -^, Xl ,X 2 ,0 [g a (xi)9a(x 2 )Pl a (0 + fa(xi)fa{x2)P la (0]}, (23) 



where l a = \n a + 1/2 1 - 1/2, l a = 2j a - l a , and g a (x) and 
f a (x) are the upper and the lower radial components of the 
reference-state wave function i/v(x)- The integration over ui 
m Eq. (|23l» can be carried out analytically in terms of the expo- 
nential integral function, as described in AppendixlAl leaving 



a 3-dimensional integration over the radial variables to be per- 
formed numerically. 



5 



III. NUMERICAL EVALUATION 

The numerical evaluation of the self-energy correction 
within the present scheme is in many respects similar to 
that in the standard potential-expansion approach. Since the 
potential-expansion method is well documented (see, e.g., a 
detailed description in Ref. I19l0 . here we concentrate on 
novel features of our evaluation as compared to the standard 
approach. They appear in the calculations of (f) the high- 
energy part of the many-potential remainder term AE™^ 
and (if) the many -potential subtraction term AE£~^ y . 

The radial integrations over x% and in the remainder 
term AE™^ y are performed after the change of variables 
(xi,x 2 ) -> \r,y) H: 



r = min(o;i, S2)/max(n, X2) , y = 2yl -£ a 2 i2- (24) 

Numerical evaluation of the radial integrals is complicated 
[specifically, for small values of Re(w)] by the presence of 
the function (E + CI) in the integrand. To explain this, 
we recall that the analytical behavior of G^°'(E + ft) is 
governed by the parameter c' = y/1 — (E + CI) 2 . Since 
E = e a — w = e a — eo — iw in the high-energy part (w G E), 
the energy argument is 



E + Cl = e a - e 



2Za 

X\ + X 2 



(25) 



For certain values of x\ and x%, Re (E + O) = 1. When 
w is small, a fast change of the phase of the square root 
-^/l — (E + Cl) 2 occurs in the vicinity of this point, which 
can lead to a numerical instability of the radial integrations. 
This problem was handled by breaking the integration inter- 
val at the point where Re (E + Cl) = 1 and employing a larger 
number of integration points in this region. 

The numerical evaluation of the subtraction term AE£^ y 
consists of a 3-dimensional integration over the radial vari- 
ables, which has a structure of the standard two-electron inte- 
gral, 



J = 



dxi dx2 



1 ^(^£2)^ ^ (26) 

-1 X12 



where the function / has a finite limit for X12 — ► 0. The inte- 
grable singularity in this expression is removed by employing 
the perimetric coordinates IB8I1 . 



U = Xl + X2 — X\2 , 
V = Xl — X 2 + X\2 , 
W = —Xl + X2 + X\2 . 



(27a) 
(27b) 
(27c) 



In the new variables, the integral J is 
1 



J = 



du dv dw x\X2 f(xi, X2, £) ■ (28) 



Performing the integrations in this expression numerically, 
one should have in mind that the function / contains a square 



root, whose argument changes its sign for certain combina- 
tions of the radial variables, similarly to the case described for 
the remainder term A.E™™ y . The point at which the argument 
of the square root vanishes is 



£0 



2Za 

Xl + x 2 



= 1 



(29) 



This feature was taken into account by breaking the integra- 
tion intervals at the singular point and by employing a larger 
number of integration points in its vicinity. 



IV. RESULTS AND DISCUSSION 

In Tables |J] [H] and [HI] we present a comparison of two 
different schemes for the evaluation of the self-energy cor- 
rection for the Is, 2s, and 2p\n states. The labels "A" and 
"B" stand for the subtraction scheme introduced in this work 
and for the standard potential-expansion approach, respec- 
tively. The entry "Free" denotes the sum of the zero- and 
one-potential terms (this part is the same in both methods), 
"Subtraction" stands for the many-potential subtraction term 
A-^many (absent in the standard approach), whereas the in- 
dividual partial-wave expansion contributions correspond to 
the many-potential remainder term AE™™y and to the many- 
potential term AE^^y in the "A" and "B" schemes, respec- 
tively. The entry "Behavior" indicates the approximate depen- 
dence of the terms of the partial-wave expansion on \k\ in the 
region of interest, i.e., for \k\ = 10 — 30. The numbers in 
parentheses represent the uncertainties in the last digit. If no 
uncertainties are indicated, numerical values are believed to 
be accurate to all digits specified. Our results obtained within 
the two approaches are compared with the numerical values 
byMohrjl. 

The comparison of the data listed in the tables demonstrates 
that the additional subtraction introduced in this work leads 
to a significant improvement of the convergence properties of 
the partial-wave expansion in all the cases studied. It also 
indicates that the new approach is applicable for the evaluation 
of the self-energy correction in the low-Z region, where the 
standard potential-expansion approach fails to yield accurate 
results. 

In the low-Z region, one has to deal with numerical can- 
cellations between individual contributions to the self-energy 
correction. The origin of these cancellations are spurious 
terms of order a(Za) 2 In Za that appear in the Feynman 
gauge when the self-energy correction is separated into the 
zero, one, and many-potential terms fiHl and that have to be 
cancelled numerically in order to obtain the physical contri- 
bution to order a(Za) 4 . In our approach, the numerical inte- 
grations can be relatively easily performed up to a sufficient 
accuracy, so that the numerical cancellations do not pose any 
serious problems. Even in the most difficult case, Z = 1, the 
present numerical scheme yields a result with a reasonable ac- 
curacy, Fi s (la) = 10.316 85(10), which is in a good agree- 
ment with the most precise value by Jentschura et al. Il ill . 
Pi, (la) = 10.316793 650(1). 
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In Table II VI we present the numerical results for the self- 
energy correction in the region that was not previously tab- 
ulated in the literature, 5 < Z < 10, and compare our nu- 
merical values for Z = 5 and 10 with evaluations by other 
authors. It is noteworthy that unlike the previous calculations 
summarized in Table IIVI our evaluation is computationally 
very cheap. The time of the calculation for one value of Z 
is less than lh on a modern personal computer. This feature 
makes the present approach very promising for extensions to 
the higher-order self-energy corrections. 

To sum up, we have developed a highly efficient scheme for 
the evaluation of the one-loop self-energy correction for an 
electron bound in a symmetric local potential (not necessarily 
the Coulomb one). The approach presented inherits the attrac- 
tive features of the standard potential-expansion method but 
yields a much better convergence rate for the resulting partial- 
wave expansion. As a result, the applicability of the potential- 
expansion method is extended into the region of large values 
of the parameter n 3 /(Za) 2 . We expect that the approach de- 
veloped will allow one to significantly improve accuracy of 
evaluations of the self-energy correction to the hyperfine split- 
ting and of the screened self-energy correction in the low-Z 
region and could be also applied for higher-order self-energy 
corrections. 



where c = y/l — (e — u>) 2 , Ch = (so — ioo, £o — i0 ] + [eo + 
iO, So + ioo), So > is the parameter of the contour, e is either 
e a or e a + f2, and the functions /j are: 

/iM = 1, (A2) 

/ 2 (w) = £12^^, (A3) 
c 

/s(w) = c, (A4) 

U(lj) = e-w, (A5) 

isH = £12^^. (A6) 

c 

Let us evaluate, e.g., the integral J\. Introducing the new 
variable y by ui = Eq + iy (ui — Eq — iy) in the upper (lower) 
part of the contour, we obtain 

Ji = — 2 Re exp[ieo £12] 

/>oo 

x dy exp{-[y + y/l + (y + ia) 2 ] xi 2 } , (A7) 
Jo 

where a — e — eq. This integral is evaluated by introducing 
the new variable 
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t = y + \/l + (y + ia) 2 - \/l-a 2 , 

with the result 

Ji = — Re exp[(ieo — a) X12] 

1 £12 exp(zx 12 ) E 1 (zxi 2 ) 

X12 z 



(A8) 



(A9) 



APPENDIX A: INTEGRALS OVER THE VIRTUAL PHOTON 
ENERGY IN THE SUBTRACTION TERM 



where a = \/l — a 2 , z = a + ia, and E± (z) is the exponential 
integral function. The results for other basic integrals are: 



The integral over ui in Eq. d23l > can be expressed as a linear 
combination of the basic integrals Ji, 



Ji = i 



dui fi(u>) exp[(i \uj\ — c) x 



12 



(Al) 



J2 = — Im exp[(ko — a) £12] 

1 — + x 2 2 exp(zxi 2 ) Ex{zx\i) 



(A10) 



J3 = — - Re exp[(ie - a) £12] 



1 Z 1 X\i I Xt, , 

YxJ 2 + ~ 2 + 2* - -z + ( 2 + -f ) eXp (^ ) E ^ 



(All) 



J 4 = - - Im exp[(ie - a) ar 12 ] 



1 Z 1 X\2 X? 2 
— + 7T~2" + o T" ex P ^12) ^l(^12 

x{ 2 %12 2z 2 z 2 



(A12) 



J5 = 



--Re exp[(ie - a) x i2 ] 



1 2 



4ia- — - r + z _fi| + ^2 +XJ3 i 2 -Z£L \ exp(za;i2) Ei(zx 12 ) 



X\2 Z 



2z 2 2z 



12 



(A13) 



All the expressions for the integrals Ji can readily be evalu- 



ated numerically. A detailed description of an algorithm for 
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TABLE I: Individual contributions to the one-loop self-energy correction for the Is state, in units of F(Za) 
new subtraction scheme, whereas "B" indicates the standard potential-expansion approach. 



"A" denotes the 





Z = 5 




Z = 10 






Z = 92 




A 


B 


A 


B 


A 


B 


Free 


-767.728 001 


767.728 


-184.021481 


184.021 48 


-0.171545 


-0.171545 


Subtraction 


30.582 424 




11.527613 




0.290 350 




\k\ = 1 


739.691981 


759.830 8 


175.775 040 


183.505 51 


1.371 144 


1.632 207 


2 


3.435 185 


8.855 9 


1.260 600 


3.339 31 


-0.001514 


0.012 042 


3 


0.227 353 


2.299 5 


0.094 384 


0.863 90 


0.001728 


0.008 313 


4 


0.029 960 


1.028 4 


0.012 972 


0.367 77 


0.000 469 


0.003 806 


5 


0.007 001 


0.568 2 


0.002 982 


0.193 45 


0.000155 


0.001988 


6 


0.002 592 


0.352 


0.001033 


0.114 57 


0.000 062 


0.001 158 


7 


0.001 246 


0.234 7 


0.000 457 


0.073 33 


0.000 029 


0.000 731 


8 


r\ nr\n coo 


0. 164 8 


0.000 231 


0.049 60 


0.000 015 


0.000 490 


9 


0.000 403 


0.120 2 


0.U00 127 


0.034 99 


0.000 008 


0.000 344 


10 


0.000 250 


0.090 4 


0.000 073 


0.025 52 


0.000 005 


0.000 251 


11 


0.000162 


0.069 7 


0.000 044 


0.019 12 


0.000 003 


0.000188 


12 


0.000108 


0.0548 


0.000 028 


0.014 65 


0.000 002 


0.000 145 


13 


0.000 074 


0.043 8 


0.000 018 


0.011 44 


0.000 001 


0.000 114 


14 


0.000 052 


0.035 5 


0.000 012 


0.009 08 


0.000 001 


0.000 091 


15 


0.000 037 


0.029 2 


0.000 008 


0.007 32 


0.000 001 


0.000 074 


S 1*1=16 


0.000 115 


0.1678 


0.000 020 


0.038 82 


0.000 002 


0.000 420 




0.000 003(2) 


0.034(3) 


0.000 001(1) 


0.007 2(4) 


0.000 000 


0.000 099(3) 


Total 


6.251627(2) 


6.252(3) 


4.654 162(1) 


4.654 1(4) 


1.490 916 


1.490 916(3) 


Ref. [6J 


6.251627(8) 




4.6541622(2) 




1.4909160(3) 


Behavior 


30/|k| 5 


100/|k| 3 


100/|k| 6 


25/|k| 3 


0.5/|k| 5 


0.25/|k| 3 



TABLE II: The same as TableQ] but for the 2s state. 







Z = 5 

A 


B 




Z = 10 

A 


B 


z 

\ 


12 


B 


Free 


-1457 


418 809 


-1457.4188 


-356 


528 846 


-356.528 8 


-1.962 337 


-1 


962 337 


Subtraction 


31 


058 101 




11 


890 558 




0.275 605 






\k\ = 1 


1410 


715 203 


1429.146 6 


339 


733 982 


346.798 4 


3.548 480 


3 


796 632 


2 


16 


099 356 


20.342 3 


6 


979 535 


8.564 2 


0.228 968 





201 262 


3 


3 


688 342 


5.475 2 


1 


722 985 


2.395 


0.073 832 





078 764 


4 


1 


379 780 


2.5818 





646 882 


1.119 1 


0.023 291 





034 477 


5 





563 792 


1.5028 





263 506 


0.643 3 


0.007 670 





017 634 


6 





233 686 


0.9791 





108 871 


0.4140 


0.002 586 





010 043 


7 





096 744 


0.685 1 





044 967 


0.286 3 


0.000 889 





006 184 


8 





039 919 


0.5039 





018 535 


0.208 1 


0.000 312 





004 042 


9 





016 454 


0.384 5 





007 645 


0.1571 


0.000114 





002 770 


10 





006 808 


0.3019 





003 174 


0.122 


0.000 044 





001 972 


11 





002 852 


0.242 5 





001 340 


0.096 9 


0.000 018 





001 449 


12 





001 228 


0.1985 





000 585 


0.078 5 


0.000 008 





001 094 


13 





000 555 


0.1650 





000 270 


0.064 5 


0.000 004 


{) 


000 845 


14 





000 271 


0.1390 





000 136 


0.053 8 


0.000 003 





000 665 


15 





000 147 


0.1184 





000 076 


0.045 3 


0.000 002 





000 533 


V 35 





000 400 


0.839 6 





000 207 


0.297 8 


0.000 005 





002 851 







000 036(10) 


0.35(8) 





000 011(7) 


0.082(10) 


0.000 000 


{) 


000 62(3) 


Total 


6 


484 865(10) 


6.54(8) 


4 


894 417(7) 


4.898(10) 


2.199 494 


2 


199 49(3) 


Ref. [6J 




6.484 8(2) 






4.89445(6) 




2.199493 8(3) 


Behavior 


5/H 4 


100/|k| 25 


10/|ft| 4 ' 5 


40/|k| 2 - 5 


io/kl 6 




2/|«| 3 



the computation of the exponential integral function of a com- plex argument can be found in Ref. Q . 
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TABLE III: The same as TableU but for the 2p 1/ / 2 state. 





Z = 5 




Z — 10 




z 


= 92 




A 


B 


A 


B 


A 


B 


Free 


-1520.728 283 


1520.728 3 


-377.853 426 


377.853 4 


-3.966 890 


-3.966 890 


Subtraction 


14.376 901 




6.031 862 




0.094 695 




|k = 1 


1481.688 696 


1483.524 8 


361.003 647 


361.809 2 


3.886 215 


3.910 643 


2 


18.863 382 


21.1676 


8.100 485 


9.121 6 


0.219 596 


0.221 754 


3 


3.816 241 


5.880 8 


1.740 929 


2.651 9 


0.058 070 


0.074 793 


4 


1.186 255 


2.832 7 


0.546 763 


1.267 1 


0.017 578 


0.032 082 


5 


0.417 922 


1.674 6 


0.193 629 


0.737 3 


0.005 943 


0.016 268 


6 


0.155 297 


1.103 9 


0.072 407 


0.477 4 


0.002 214 


0.009 235 


7 


0.U59 4 f 1 


0.779 4 




o.oou y 


n ruin mo 


0.000 686 


,s 


0.023 323 


0.577 2 


0.011 192 


0.240 5 


0.000 420 


0.003 724 


9 


0.009 401 


0.442 7 


0.004 654 


0.181 1 


0.000 215 


0.002 560 


10 


0.003 945 


0.348 9 


0.002 053 


0.140 2 


0.000 122 


0.001830 


11 


0.001764 


0.281 1 


0.000 985 


0.1109 


0.000 074 


0.001350 


12 


0.000 866 


0.230 4 


0.000 528 


0.089 3 


0.000 048 


0.001 024 


13 


0.000 479 


0.191 7 


0.000 318 


0.073 


0.000 033 


0.000 794 


14 


0.000 300 


0.161 5 


0.000 212 


0.060 5 


0.000 023 


0.000 628 


15 


0.000 209 


0.137 5 


0.000 153 


0.050 7 


0.000 016 


0.000 505 


£|*|=le 


0.000 928 


0.9610 


0.000 668 


0.320 3 


0.000 054 


0.002 747 


£ M=36 ( eXtr -) 


0.00013(3) 


0.32(8) 


0.000 08(2) 


0.079(10) 


0.000 002 


0.000 61(3) 


Total 


-0.122 77(3) 


-0.11(8) 


-0.11484(2) 


-0.113(10) 


0.319 341 


0.319 34(3) 


Ref. [6J 


-0.1228(2) 




-0.11483(4) 




0.319 340 8(4) 


Behavior 


2/M 3 - 5 


130/|k| 25 


7/|k| 4 


200/|k| 3 


10/|k| 5 


2/|k| 3 



TABLE IV: One-loop self-energy correction, in terms of F(Za). 





Is 


2s 


2pi/ 2 


2P3/2 


Ref. 
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0.125 6(1) 
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[11] 
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